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^ : Abstract 



Thermal effects on the creation of particles under the influence of time- 



O^ . dependent boundary conditions are investigated. The dominant temperature 
correction to the energy radiated by a moving mirror is derived by means of 

^ \ response theory. For a resonantly vibrating cavity the thermal effect on the 

(N . 

Psj . number of created photons is obtained non-perturbatively. Finite tempera- 

\^ ' tures can enhance the pure vacuum effect by several orders of magnitude. The 

Qs \ relevance of finite temperature effects for the experimental verification of the 



dynamical Casimir effect is addressed. 



Q^; PACS-numbers: 42.50.Lc; ll.lO.Wx; ll.lO.Ef; 03.70.+k 

-(— > 

j^ ■ In 1948 H. B. G. Casimir predicted an attractive force between two perfectly conducting, 

qh| parallel plates placed in the vacuum [0. A variety of fundamental and measurable conse- 

quences of quantum fluctuations under the influence of external conditions have been derived 
/\ • during the last decades (see e.g. 0] for reviews). With the recent precision measurement of 

. S : the Casimir force performed by Mohideen et at B and by Lamoureaux § conclusive test, 

of Casimir's prediction are now available. They confirm the basic concepts of quantum field 
theory in the presence of static external constraints. 

One could expect that the success in measuring the static Casimir force could also in- 
tensify experimental efforts in order to verify a not less fundamental prediction, namely 
the dynamical Casimir effect, i.e., the creation of particles out of the vacuum induced by 
the interaction with dynamical external constraints. The phenomenon of vacuum radiation 
induced by moving mirrors [Q initiated intensive studies (see e.g. and references therein). 
In particular the creation of photons in vibrating cavities seems to be the most promising 
scenario for a possible experimental verification of motion-induced vacuum radiation (see 

e.g. wm- 



The thermodynamics of the static Casimir effect has been investigated intensively (see 
e.g. |13[ and references therein). In that context temperature effects are known to even 



dominate the pure vacuum effect (at T = 0) and, in consequence, have to be taken into 
account when analysing the data in measurements of the static Casimir force. In contrast 
to this the dynamical Casimir effect at finite temperature so far has not been subject of 
research. It has been anticipated in a recent investigation by Lambrecht et al. [Q that 
temperature effects could play an important role for the generation of a photon pulse in a 
vibrating Fabry-Perot cavity. However, realistic calculations of thermal effects on quantum 
radiation within the framework of quantum field theory of time-dependent systems at finite 
temperature are not yet available. 

Accordingly, it is our major intention to close this gap and to provide a generalization 
of the Hamiltonian approach presented recently in Refs. 0,^]. In this letter we focus 
the discussion predominantly on results obtained for the thermal contribution to photon 
production in a resonantly vibrating cavity as one of the most relevant configuration when 
aiming for experimental tests of the dynamical Casimir effect [Q. We like to address the 
question whether or not the effect of quantum radiation might be covered by the thermal 
background and we will examine the conditions under which it remains most significant even 
at large temperatures. Details of the derivations and further applications of the formalism 
will be presented in a forthcoming publication [ p^ . 

We recall the canonical formalism in Ref. []n[, where we considered a constrained, non- 



interacting, real, massless scalar field. Similarly, the formalism also holds for bosonic quan- 
tum fields interacting with classical external background fields [|15|. The boundary respec- 
tively the background may undergo small but arbitrary dynamical changes resulting in an 
additional interaction Hamiltonian Hi which is assumed to be switched on and off at asymp- 
totic times t -^ — oo and t -^ +00, respectively. For the boundary being initially at rest the 
(closed) system consisting of the scalar field enclosed by the boundary is assumed to be at 
thermal equilibrium described by a statistical operator /)(t -^ —00) = p^. The existence of 
a total Hamiltonian H = Hq + Hi describing the evolution of the system implies that any 
kind of backreaction of the quantum field upon the dynamics of the external constraints or 
any relaxation processes will be neglected. In addition, there should be no measurements on 
the quantum systems during the dynamical phase. We adopt the interaction representation. 
When the boundary or the background field experiences dynamical changes the system will 
no longer remain at thermal equilibrium. The time-evolution of the statistical operator p is 
governed by the time-dependent interaction Hamiltonian Hi resulting in the von Neumann 
equation (units where h = c = ks = I are used throughout): 



^^, = [Hi,p\ , (1) 

together with the initial condition p{t —>■ — oo) = p^. Here we disregard an exphcit time- 
dependence {dp/dt)exp. that could account for relaxation or backreaction processes. Accord- 
ingly, this equation of motion can be integrated formally with the aid of the time-evolution 
operator 



U{t',t) =T 



exp —i / dti Hi{ti 



(2) 



where T denotes time-ordering. The thermal expectation value at asymptotic times t -^ oo 
of any relevant observable A is determined by 

(i) = Tr{ip(i^oo)} = Tr{if/pof/+} ■ (3) 

The trace Tr{- ■ ■} involved is taken most conveniently over the Fock space of the unper- 
turbed Hamiltonian Hq refering to the initial ensemble (with Tr{po} = 1). The microscopic 
entropy S = — Trjplnp} = — Trjpolnpo} = 'S'o remains constant in time. A more detailed 
examination of the thermodynamical aspects involved (e.g. effective entropy) will be given 
T6| . Eqs. (|I]) - d^) generalize the canonical approach to quantum radiation of Ref . p!^J15[ 



m 



to account for finite-temperature effects as well. In order to explore the influence of finite 
temperatures on the dynamical Casimir effect we have to investigate the thermal expecta- 
tion value of the number operator Nx = al^ax of particles with frequency Q^. Note, that the 
proper definition of (quasi-) particle creation- (annihilation-) operators a^ {d\) is provided 
with respect to the ground state |0) of the unperturbed Hamiltonian Hq, which becomes 
diagonal, i.e., Hq = Ql {N\ + 1/2). The total radiated energy associated with the dynamical 
Casimir-effect can be deduced from the expectation value Eq. (|^) of the normal-ordered 
unperturbed Hamiltonian Tt{:Hq: p}. 

In the following we consider the generic case of a constrained massless scalar or vector 
field at finite temperature T = 1//? described initially (t -^ — oo) by the statistical operator 
(canonical ensemble) 



Po 



exp(-(3Ho)/Zo . (4) 



The closed system is assumed to be initially at thermal equilibrium. The system leaves 
the thermal equilibrium when undergoing some dynamical changes. Eq. (^ allows for 
a systematic perturbative approach in those cases, for which a closed expression for the 
time-evolution operator is not available. Especially for the thermal expectation value of the 



spectral number density {Nx) one obtains up to quadratic response neglecting corrections 
of order 0{Hf) 

(iVA) = Tr{iVAPo} + Tr|iVA y dt Hi{t) , po] J dt H,{t)^ 

= {Nx)o + ANx . (5) 

The first term {Nx)o = Tt{NxPo} represents the usual thermal Bose-distribution function 
while the second term AA^^ denotes the motion-induced change of the number of particles 
(e.g. photons) with frequency fix ^^ ^ given temperature T. The total radiated energy, 
i.e. the change of the internal energy associated with quantum radiation is determined by 
AE = ^lANx and can be deduced according to E = {-.Hq:) = (:^o:)o + ^E. 

Let us now turn to the generic particle-creation processes. In the case of constrained 
bosonic quantum fields satisfying linear equations of motion the general form of the self- 
adjoint perturbation Hamiltonian induced by the interaction with the external conditions 
assumes the rather general form (sum convention) 

r ^ 1 

J dt Hi{t) = - (S^^ala^ + Sl^d^^d^^ + U^^al^dy + C (6) 

with S^y = Sy^ and lA^y = W*^,. The constant C gives rise to a pure phase factor and thus 
drops out in any expectation value (^. The iS-term may be interpreted as a generator of 
a multi-mode squeezing operator and the W-term may be envisaged as a hopping operator. 
Evaluating the quadratic response according to Eq. (^ the change of the number of particles 
yields 

ANx = \Sx/{l + {N,)o + {Nx)o) 

+ \Ux/ ((iV,)o - {Nx)o) . (7) 

Note, that the iS-term contains the pure vacuum contribution AA^A(vac) = |iSap| . One 
observes that only the iS-term contributes to the total number of created particles (sum over 
A), while the W-term does not increase the total number of particles (similar terms appear 
in master-equations). However, it modifies the configurations of occupied particle states 
within the ensemble and thus also increases the total energy. 



For trembling cavities the perturbation Hamiltonian (cf. Ref. [14| for details) 



H,it) = qlit) AfiJ(t)/2 + q,p,M,,{t) (8) 

appears as a sum of a squeezing (first term) and a velocity term (second term). The squeezing 
term involves the deviations AVti of the eigenmode spectrum induced by a change of the 



shape of the cavity, while the velocity contribution involves the (antisymmetric) intermode 
couplings A4ij,u arising from the motion of the boundary. 

In scenarios, where the velocity effect is supposed to be negligible the matrix S be- 
comes diagonal and the W-term does not contribute at all. In view of Eq. (^ we 
identify the corresponding particle-production rate as the product of the pure vacuum- 
squeezing effect AA^f (vac) = \Sxx\ times a thermal distribution factor, i.e. AA^f = 
AiV^(vac) (l + 2(iVA)o)- It gives rise to a linear dependence on T in the high-temperature 
limit. 

The pure velocity effect may be illustrated most simply by considering a single moving 
mirror in 1+1 dimensions. For a single mirror placed at a position 77 (t) the discrete eigen- 
modes fl^ have to be replaced by the continuous variable k. The total radiated energy is 
derived as 

E = -^jdtf{t) + ^-T' jdtfi'it) , (9) 

which generalizes the zero-temperature result (first term) obtained by Fulling and Davis 
and by Ford and Vilenkin 0. The ratio of the finite temperature correction to the radiated 
energy and the pure vacuum contribution turns out to be of the order (9(T^r^), where r 
denotes a characteristic time scale of the underlying dynamics. 

Let us now investigate the finite-temperature effects on the dynamical Casimir-effect in 
a resonantly vibrating cavity. In order to allow for an experimental verification the number 
of motion-induced created particles should be as large as possible. One way to achieve 
this goal is to utilize the phenomenon of parametric resonance, which occurs in the case of 
harmonically time-dependent perturbations characterized by some frequency u. Obtaining 
large numbers may indicate that one has left the region, where second-order perturbation 
theory does apply. 

In the case of oscillating disturbances, however, it is possible to evaluate the time- 
evolution operator to all orders of Hi analytically employing yet another approximation, 
the so-called rotating wave approximation (see e.g. [^). Let us assume that the explicit 
time-dependence of the perturbation Hamiltonian possesses an oscillatory behaviour like 
e sin(2c(jt) during a sufficiently long period of time T, such that the conditions uoT ^ 1, 
£ <^ 1 and eujT = 0{1) hold. The main consequence of the rotating wave approxima- 
tion consists in the fact that it allows for the derivation of a time-independent, effective 
Hamiltonian H^ after performing the integration over time according to Eq. (^: 

dtHi^THf . (10) 





This effective Hamiltonian accounts essentially for the dominant contribution arising from 
the frequency uj. For the evolution operator Eq. (^ this approximation implies the neglec- 
tion of all time-ordering effects (all oscillating terms involving commutators average out). 

For a vibrating cavity the effective Hamiltonian H^ is easily calculated from the inter- 
action operator (^. Assuming a harmonic time-dependence proportional to esm{2ujt) or 
ecos{2ujt) for both, the squeezing (Afi^) and the velocity terms {M.^y) only those terms will 
survive, which match the resonance conditions. For the squeezing term the resonance condi- 
tion reads Vl\ = uj and for the velocity term \Vt^ ± i7[J| = 2u;, respectively. Accordingly, the 
squeezing effect creates always particles with the frequency u provided this cavity mode does 
exist. We restrict our further consideration to the situation, where the oscillation frequency 
UJ corresponds to the lowest cavity mode uj = min {fi^} = fi?. The fundamental resonance 
frequency is determined by the characteristic size L of the cavity, e.g. fij = v^vr/L for a 
cubic cavity. For the lowest mode the resonance condition for squeezing Q^ = Q^ is satisfied 
automatically. Whether the resonance condition for the velocity effect |fij^ ± r2°| = 2Q^ can 
be satisfied or not depends on the spectrum of the particular cavity under consideration. For 
a one- dimensional cavity the eigenvalues Q^ are proportional to integers and thus it can be 
satisfied leading to an additional velocity contribution. For most cases of higher-dimensional 
cavities, e.g. a cubic one, this condition cannot be fulfilled. Thus the velocity effect does 
not contribute within the rotating wave approximation (cf. Ref. 0). 

In consequence, only the squeezing term contributes. The effective Hamiltonian can be 
derived immediately for the lowest cavity mode. From the contributing Afi^-terms one 
obtains Hf^ = iD^le/A (aj)^ — (tii)^ respectively for the time-evolution operator 

f7(T,0)^exp|^[(al)2-(aO^]|=5i . (11) 

Note, that within this approximation U coincides with a squeezing operator Si for the lowest 
mode A = 1. This confirms the notion of the Ai7^-terms in (^ as squeezing contribution. 
After having derived a closed expression for the time-evolution operator this enables us to 
calculate the expectation value for the number operator to all orders in Hf^ 

iV.)-(iV.)^ + 5.isinh^(£^)(l + 2(iV,)J . (12) 

This non-perturbative result implies that at finite temperature the number of photons AA^i 
created resonantly in the lowest cavity mode increases exponentially. The vacuum creation 
rate AA'^f (vac) = sinh (erilT/2) (see Ref. [§) gets enhanced by a thermal distribution factor. 



Equation (|T2D represents one essential result of our investigations. In order to indicate 
its experimental relevance one may specify the characteristic parameters. In Fig. |^ the total 
number of photons in the fundamental mode u = Q^ inside of a cavity of a characteristic 
size L is depicted as a function of the time-duration T of the vibration of its wall. The 
contribution AA^^i of additional, motion-induced photon creation at T = 0, i.e. the pure 
vacuum effect (dashed curve) is compared with the corresponding contribution at a temper- 
ature T (solid curve). The shadowed area indicates the number of photons (iVi)o present 
at the temperature T before the system will undergo dynamical changes together with its 



thermal variance ao{Nx) = V (iV|)o - (A^a )o ~ (^a )o for {Nx )o > 1- The latter reflects the 
uncertainty when measuring the number of photons at a given temperature T. For room 
temperature T ^ 290 K, which corresponds to thermal wave lengths of about 50 fim and 
considering a cavity of a typical size L ^ 1 cm one obtains a thermal factor (1 + 2(A'^i )o) of 
order 10^. As a consequence after the vibration time T the number of photons AA^^i created 
by the dynamical Casimir effect at the given temperature will be three orders of magnitude 
larger in comparison with the pure vacuum effect at T = 0. At finite temperatures the 
dynamical Casimir effect should become observable even after shorter vibration times T as 
one would expect from looking at the pure vacuum effect. This strong enhancement of the 
dynamical Casimir effect occuring at finite temperatures could be exploited in experiments 
to verify the phenomenon of quantum radiation as long as backreaction processes can be 
neglected. Of course, one has to ensure conditions that will lead to a significant vacuum 
effect as well. This implies that the argument of the hyperbolic sine function in Eq. (0), i.e. 
the squeezing parameter eujT/2 should be at least of order 1. An estimate of the maximal 
value of the dimensionless amplitude of the resonance wall vibration emax < 10~^ is given 
in Ref. 0. It still remains as a challenge whether or not the requirement £:ti;T/2 ^ 1 could 
be achieved in a realistic experiment. We conclude: Provided an experimental device for 
generating a considerable vacuum contribution becomes feasible, then we predict a strong 
enhancement of the dynamical Casimir effect at finite temperature. From the theoretical 
point of view there is no definite need to perform an even more involved experiment at low 
temperatures. 
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FIG. 1. Number of photons A'^i produced via the dynamical Casimir effect inside a resonantly 
vibrating cavity for a fixed value of the squeezing parameter eujjl as a function of the time duration 
T. The fundamental frequency is choosen to be w = 146 GHz corresponding to a typical size of 
the cavity of 1 cm. For the dimensionless amplitude a value of e = 6 • 10"^*^ has been assumed. A 
finite temperature (T = 290 K) can enhance the effect (solid line) by several orders of magnitudes 
compared with the pure vacuum effect at T = (dashed line) . The thermal background is indicated 
by the shadowed area (half of the thermal variance uq of the photon number A^i). 



